Universal behavior of extreme value statistics for selected observables of dynamical 

systems 
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The main results of the extreme value theory developed for the investigation of the observables 
of dynamical systems rely, up to now, on the Gnedenko approach. In this framework, extremes 
are basically identified with the block maxima of the time series of the chosen observable, in the 
limit of infinitely long blocks. It has been proved that, assuming suitable mixing conditions for 
the underlying dynamical systems, the extremes of a specific class of observables are distributed 
according to the so called Generalized Extreme Value (GEV) distribution. Direct calculations 
show that in the case of quasi-periodic dynamics the block maxima are not distributed according 
to the GEV distribution. In this paper we show that, in order to obtain a universal behaviour 
of the extremes, the requirement of a mixing dynamics can be relaxed if the Pareto approach is 
used, based upon considering the exceedances over a given threshold. Requiring that the invariant 
measure locally scales with a well defined exponent - the local dimension -, we show that the limiting 
distribution for the exceedances of the observables previously studied with the Gnedenko approach 
is a Generalized Pareto distribution where the parameters depends only on the local dimensions 
and the value of the threshold. This result allows to extend the extreme value theory for dynamical 
systems to the case of regular motions. We also provide connections with the results obtained with 
the Gnedenko approach. In order to provide further support to our findings, we present the results 
of numerical experiments carried out considering the well-known Chirikov standard map. 



I. INTRODUCTION 

Extreme value Theory was originally introduced by 
Fisher and Tippett [1] and formalised by Gnedenko [2], 
who showed that the distribution of the maxima of a sam- 
ple of independent identically distributed (i.i.d) stochas- 
tic variables converges under very general conditions to 
a member of the so-called Generalised Extreme Value 
(GEV) distribution. The attention of the scientific com- 
munity to the problem of understanding extreme values 
theory is growing, also because this theory is crucial in a 
wide class of applications for defining risk factors such as 
those related to instabilities in the financial markets and 
to natural hazards related to seismic, climatic and hydro- 
logical extreme events. Even if the probability of extreme 
events decreases with their magnitude, the damage that 
they may bring increases rapidly with the magnitude as 
does the cost of protection against them. From a theo- 
retical point of view, extreme values of observables are 
related to large fluctuations of the corresponding under- 
lying system. An extensive account of recent results and 
relevant applications is given in [3]. 

The traditional ('Gnedenko') approach for the statisti- 
cal inference of extremes is related to the original results 
by Gnedenko [2]: we partition the experimental time se- 
ries into bins of fixed length, we extract the maximum of 
each bin, and fit the selected data to the GEV distribu- 
tion family using methods such as maximum likelihood 
estimation (MLE) or L-moments. See [4] for a detailed 
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account of this methodology. The selection of just one 
maximum in a fixed period may lead to the loss of rele- 
vant information on the large fluctuations of the system, 
especially when there are many large values in a given 
period [5]. This problem can be taken care of by con- 
sidering several of the largest order statistics instead of 
just the largest one. For such maxima distributions we 
expect convergence to the Generalized Pareto Distribu- 
tion (GPD) introduced by Pickands III [6] and Balkema 
and De Haan [7] to model the exceedances over a given 
threshold. We call this the 'Pareto' approach. Also this 
approach has been widely adopted for studying empiri- 
cally natural extreme phenomena such as those related 
to waves, winds, temperatures, earthquakes and floods 
[8-10]. 

Both the Gnedenko and Pareto approaches were orig- 
inally designed to study extreme values for series of i.i.d 
variables. In this case it is well known that a strong 
connections exists between the two methodologies, as we 
have that if block maxima obey the GEV distribution, 
then exceedances over some high threshold will have an 
associated GPD. Moreover, the shape parameter of the 
GPD and that of the corresponding GEV distribution are 
identical [11]. As a result, several practical methods (e.g. 
Hill's and Pickands' estimators) developed for estimat- 
ing the shape parameter of the GEV distribution of the 
extremes of a given time series are actually based upon 
comparing the GPD fits at various thresholds [12, 13]. In 
practical terms, it appears that, while the Gnedenko and 
Pareto approaches provide equivalent information in the 
asymptotic limit of infinitely long time series, the GPD 
statistics is more robust when realistic, finite time series 
are considered (see, e.g., [14]). 

In recent years, especially under the influence of the 
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rapid development of numerical modelling in the geo- 
physical sciences and of its applications for the investiga- 
tion of the socio-economic impacts of extreme events, it 
has become of great relevance to understand whether it 
is possible to apply the extreme value theory on the time 
series of observables of deterministic dynamical systems. 
Carefully devised numerical experiments on climate mod- 
els of various degrees of complexity have shown that the 
speed of convergence (if any) of the statistical properties 
of the extremes definitely depends on the chosen climatic 
variable of interest [4, 15-17]. 

Several papers have addressed this issue at a more gen- 
eral level. A first important result is that when a dy- 
namical system has a regular (periodic of quasi-periodic) 
behaviour, we do not expect, in general, to find conver- 
gence to GEV distributions for the extremes of any ob- 
servable. These results have been presented by Balakr- 
ishnan et al. [18], and more recently, by Nicolis et al. [19] 
and by Haiman [20] . 

A different mathematical approach to extreme value 
theory in dynamical systems was proposed in the land- 
mark paper by Collet [21], which has paved the way for 
the recent results obtained in the last few years [22-25]. 
The starting point of all of these investigations has been 
to associate to the stationary stochastic process given by 
the dynamical system, a new stationary independent se- 
quence which obeys one of the classical three extreme 
value laws introduced by Gncdenko [2]. The assump- 
tions which are necessary to observe a GEV distribu- 
tion in dynamical systems rely on the choice of suitable 
observables (specific functions of the distance between 
the orbit and the initial condition, chosen to be on the 
attractor) and the fulfillment of particular mixing con- 
ditions that guarantee the independence of subsequent 
maxima. Recent studies have shown that the resulting 
parameters of the GEV distributions can be expressed as 
simple functions of the local (around the initial condi- 
tion) dimension of the attractor, and detailed numerical 
investigations have clarified the conditions under which 
convergence to the theoretical GEV distributions can be 
satisfactorily achieved when considering finite time series 
[26-28]. 

In this paper, we wish to attempt a unification of these 
two lines of work by using the Pareto rather than the 
Gncdenko approach. We choose the same class of ob- 
servables presented in [22-28] and show that, assuming 
only that the local measures scales with the local dimen- 
sion [29], it is possible to obtain by direct integration 
a GPD for the threshold exceedences when considering 
a generic orbit of a dynamical systems, without requir- 
ing any special mixing properties. The parameters will 
depend only on the choice of the threshold and, more 
importantly, on the local dimension. Note that Castillo 
and Hadi [5] had already pointed out that in the case 
of periodic or quasi-periodic motion the Gnedenko ap- 
proach to the evaluation of the extreme value statistics 
is inefficient, basically because in the limit of very large 
blocks, we tend to observe always the same maximum 



in all bins. To support our analytical results we provide 
numerical experiments that we carry out considering the 
classic Chirikov standard map [30] . This paper is organ- 
ised as follows. In section 2 we recapitulate the extreme 
value theory for dynamical systems obtained using the 
Gncdenko approach. In section 3 we present our general 
results obtained using the Pareto approach. In section 
4 we provide support to our investigation by examining 
the results of the numerical simulations performed on the 
standard map. In section 5 we present our final remarks 
and future scientific perspectives. 



II. GNEDENKO APPROACH: GENERALIZED 
EXTREME VALUE DISTRIBUTIONS IN 
DYNAMICAL SYSTEMS 

Gnedenko [2] studied the convergence of maxima of 
i.i.d. variables 

Xq, X\, ..X m -i 

with cumulative distribution function (cdf) F(x) — 
P{a m (M m - b m ) < x} where a m and b m are normaliz- 
ing sequences and M m = max{V , Xi, X m -i}. Under 
general hypothesis on the nature of the parent distribu- 
tion of data, Gncdenko [2] showed that the asymptotic 
distribution of maxima, up to an affine change of vari- 
able, belongs to a single family of generalized distribu- 
tion called GEV distribution whose cdf can be written 
as: 



F GEV (x;^,a,K) = c 



-t(x) 



(1) 



where 



= J(1 + K (^))- 1/K if^O 

e -(x- M )/a if K = 



This expression holds for 1 + k(x — n)/a > 0, using 
(location parameter) and a > (scale parameter) as scal- 
ing constants in place of b m , and a m [31], in particular, 
in Faranda et al. [27] wc lia,v6 shown that fi — b rn and 
a = l/a m , where k g 1 is the shape parameter (also 
called the tail index). When the distribution cor- 

responds to a Gumbel type (Type 1 distribution). When 
the index is positive, it corresponds to a Frechet (Type 2 
distribution); when the index is negative, it corresponds 
to a Weibull (Type 3 distribution). 

Let us consider a dynamical system (f2, B, v, /), where 
fl is the invariant set in some manifold, usually M. d , B 
is the Borel a-algebra, / : il — > Vl is a measurable map 
and v an /-invariant Borel measure. In order to adapt 
the extreme value theory to dynamical systems, follow- 
ing [22-25], we consider the stationary stochastic process 
X ,X 1} ... given by: 



X m (x)=g(dist(f m (x),()) VmeN 



(3) 
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where 'dist' is a distance in the ambient space f2, £ is a 
given point and g is an observable function. The partial 
maximum in the Gncdcnko approach is defined as: 



53-type observable: 



M„ 



max{X , ...,X TO _i}. 



(4) 



Defining r = dist(x, C), we consider the three classes of 
observables gi, i = 1,2,3: 



(5) 
(6) 

(7) 



gi (r) = -log(r) 
g 2 (r) = r~ p 
g 3 (r) = C-rP 

where C is a constant and (3 > £ R. Using the ob- 
servable gi we obtain convergence of the statistics of the 
block maxima of their time series obtained by evolving 
the dynamical system to the Type i distribution if one 
can prove two sufficient conditions called D 2 and D', 
which basically imply a sort of independence of the series 
of extremes resulting from the mixing of the underlying 
dynamics [22]. The conditions cannot be simply related 
to the usual concepts of strong or weak mixing, but are 
indeed not obeyed by observables of systems featuring a 
regular dynamics or power-law decay of correlations. 

A connection also exists between the existence of ex- 
treme value laws and the statistics of first return and 
hitting times, which provide information on how fast the 
point starting from the initial condition £ comes back 
to a neighborhood of £, as shown by Freitas et al. [23] 
and Freitas et al. [32]. In particular, they proved that 
for dynamical systems possessing an invariant measure 
v, the existence of an exponential hitting time statistics 
on balls around ^-almost any point ( implies the exis- 
tence of extreme value laws for one of the observables 
of type gi,i — 1,2,3 described above. The converse is 
also true, namely if we have an extreme value law which 
applies to the observables of type gi,i = 1,2,3 achiev- 
ing a maximum at £, then we have exponential hitting 
time statistics to balls with center £. Recently these re- 
sults have been generalized to local returns around balls 
centered at periodic points [24]. 

In Faranda et al. [26, 27, 28] we analised both from an 
analytical and numerical point of view the extreme value 
distribution in a wide class of low dimensional maps. We 
divided the time series of length k of the gi observables 
into n bins each containing the same number m of ob- 
servations, and selected the maximum (or the minimum) 
value in each of them [33]. We showed that at lead- 
ing order (the formulas are asymptotically correct for 
m, k — > 00), the GEV parameters in mixing maps can 
be written in terms of m (or equivalently n) and the lo- 
cal dimension of the attractor D. We have: 



(/i-type observable: 



1 1 
a = — fj, ~ — ln(fc/n) 



K = 



52-type observable: 



a ~ n D (i ~ n D K = — 



(8) 



(9) 



a ~ n D 



H = C k= - — 



(10) 



Moreover, we clearly showed that other kind of distribu- 
tions not belonging to the GEV family are observed for 
quasi-periodic and periodic motions. 



III. THE PARETO APPROACH: 
GENERALIZED PARETO DISTRIBUTIONS IN 
DYNAMICAL SYSTEMS 

We define an exceedance as z = X — T, which mea- 
sures by how much X exceeds the threshold T. As dis- 
cussed above, under the same conditions under which 
the block maxima of the i.i.d. stochastic variables X 
obey the GEV statistics, the exceedances z are asymp- 
totically distributed according to the Generalised Pareto 
Distribution [11]: 



F GPD (z:£,<7) = 



for £ ^ 0, 



l-exp(-f) for£ = 0, 



(11) 



where the range ofzis0<z<ooif£ <0 and 
0<2<cr/£if£>0. We consider the same set up 
described in the previous section and take into account 
an observables g = g(dist(x,Q) — g{r), such that g 
achieves a maximum g max for r = (finite or infinite) and 
is monotonically decreasing. We study the exceedance 
above a threshold T defined as T — g(r*). We obtain 
an exceedence every time the distance between the or- 
bit of the dynamical system and ( is smaller than r*. 
Therefore, we define the exceedances z — g(r) — T. By 
the Bayes' theorem, we have that P(r < g^ 1 {z + T)\r < 
g~ x (T)) = P(r < g- 1 (z + T))/P(r < g _1 (T)). In terms 
of invariant measure of the system, we have that the 
probability H 3 ^t(z) of observing an exceedance of at least 
z given that an exceedence occurs is given by: 



Hg,T{z) 



v{B g - 



(0) 



v(B g -i {T )(Q) 



(12) 



Obviously, the value of the previous expression is 1 if 
z = 0. In agreement with the conditions given on g, 
the expression contained in Eq. (12) monotonically de- 
creases with z and vanishes when the radius is given by 
9 \9max)- Note that the corresponding cdf is given by 
Fg.T(z) = l-Hg.T(z). In order to address the problem of 
extremes, we have to consider small radii. At this regard 
we will invoke, and assume, the existence of the following 
limit 

lim lQ g^( g ' (0) = d ,q for ^ chosen v _ / 13 ) 
r->0 logr 

where D(Q is the local dimension of the attractor [29]. 
Therefore, we rewrite the following expression for the tail 
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probability of exceedance: 
H b ,t{z) - 



D 



(14) 



where we have dropped the £ dependence of D to simplify 
the notation. By substituting g with specific observable 
we are considering, we obtain explicitly the correspond- 
ing extreme value distribution law. 

By choosing an observable of the form given by either 
gi, g 2 , or g 3 , we derive as extreme value distribution 
law one member of the Generalised Pareto Distribution 
family given in Eq. (11). Results are detailed below: 



• Pi-type observable: 



D 

• 52-type observable: 
T0 
~~D 



a = 



• <?3-type observable: 

(C - T)(3 



a = 



D 



I 
D ' 



(15) 



(16) 



(17) 



The previous expressions show that there is a simple al- 
gebraic link between the parameters of the GPD and the 
local dimension of the attractor around the point £. This 
implies that the statistics of extremes provides us with a 
new algorithmic tool for estimating the local fine struc- 
ture of the attractor. These results show that it is possi- 
ble to derive general properties for the extreme values of 
the observables g\, g 2 , or g% independently on the qualita- 
tive properties of the underlying dynamics, be the system 
periodic, quasi-periodic, or chaotic. Therefore, by taking 
the Pareto instead of the Gnedenko approach, we are able 
to overcome the mixing conditions (or the requirements 
on the properties of the hitting time statistics) needed to 
derive a general extreme value theory for dynamical sys- 
tems, as proposed in [21-25]. In [26-28] we had proposed 
that the reason why a link between the extreme value the- 
ory and the local properties of the invariant measure in 
the vicinity of the point £ can be explained by the fact 
that selecting the extremes of the observables g\ , g 2l or g^ 
amounts to performing a zoom around £. In the case of 
the Gnedenko approach, such a picture is accurate only 
if the dynamics is mixing (time and spatial selection cri- 
teria are equivalent). Instead, in the case of the Pareto 
approach, this is literally what we are doing when writing 
Eq. (14), as we are remapping the radius of the ball in a 
monotonic fashion though the inverse of the g-functions. 



A. Relationship between the Gnedenko and Pareto 
approaches 

The relation between GEV and GPD parameters have 
been already discussed in literature in case of i.i.d vari- 
ables [13, 14, 34, 35]. Coles [13] and Katz et al. 



[34] have proven that the cdf of the GEV defined as 
Fqev {z\ M; a i K ) can be asymptotically written as that 
of GPD under a high enough threshold as follows: 



F G ev(z; n,a,K) ~ F GPD (z; T, a, £) = 

= i-[i-e(^)] 1/? 



(18) 



where k = f, o = a. + £(T - fi), and T = fj,+ |(A"« - 1), 
with ln(a) = ln(er) + £ ln(A). In the present case, we have 
to compare Eqs. (8)-(10) for GEV with Eqs. (15)-(17) 
for GPD, keeping in mind that the GEV results hold only 
under the mixing conditions discussed before. While it 
is immediate to check that n = £, the other relationships 
are valid in the limit of large n, as expected. 



IV. NUMERICAL INVESTIGATION 

The standard map [36] is an area-preserving chaotic 
map defined on the bidimensional torus, and it is one of 
the most widely-studied examples of dynamical chaos in 
physics. The corresponding mechanical system is usually 
called a kicked rotator. It is defined as: 




Vt- £ sin(27ra; t ) 
■xt+yt + l 



mod 1 
mod 1 



(19) 



The dynamics of the map given in Eq. (19) can be regular 
or chaotic. For K « 1 the motion follows quasi periodic 
orbits for all initial conditions, whereas ii K >> 1 the 
motion turns to be chaotic and irregular. An interesting 
behavior is achieved when K ~ 1: in this case we have 
coexistence of regular and chaotic motions depending on 
the chosen initial conditions [37]. 

We perform for various values of K ranging from 
K = 1CP 4 up to K = 10 2 an ensemble of 200 simula- 
tions, each characterised by a different initial condition £ 
randomly taken on the bidimensional torus, and we com- 



pute for each orbit the observables gi 



1, 2, 3. In each 



case, the map is iterated until obtaining a statistics con- 
sisting 10 4 exceedances, where the threshold T = 7 ■ 10~ 3 
and (3 = 3. We have carefully checked that all the re- 
sults are indeed robust with respect to the choice of the 
threshold and of the value of 0. For each orbit, we fit 
the statistics of the 10 4 exceedances values of the observ- 
ables to a GPD distribution, using a MLE estimation 
[5] implemented in the MATLAB© function gpdfit [38]. 
The results are shown in Fig. 1 for the inferred values 
of £ and a and should be compared with Eqs. (15)-(17). 
When K <gC 1, we obtain that the estimates of £ and a 
are compatible with a dimension D = 1 for all the initial 
conditions: we have that the ensemble spread is negligi- 
ble. Similarly, for K ^> 1, the estimates for £ and a agree 
remarkably well with having a local dimension D = 2 for 
all the initial conditions. In the transition regime, which 
occurs for K ~ 1, the ensemble spread is much higher, 
because the scaling properties of the measure is differ- 
ent among the various initial conditions. As expected, 
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the ensemble averages of the parameters change mono- 
tonically from the value pertaining to the regular regime 
to that pertaining to the chaotic regime with increasing 
values of K. Basically, this measures the fact that the 
so-called regular islands shrink with K. Note that in the 
case of the observable gi, the estimate of the £ is ro- 
bust in all regimes, even if, as expected, in the transition 
between low and high values of K the ensemble spread 
is larger. These results can also be compared with the 
analysis presented in Faranda et al. [26], where we used 
the Gnedenko approach. In that case, the values ob- 
tained in the regular regions were inconsistent with the 
GEV findings, the very reason being that the dynamics 
was indeed not mixing. Here, it is clear that the statis- 
tics can be computed in all cases, and we have a powerful 
method for discriminating regular from chaotic behaviors 
through the analysis of the inferred local dimension. 



V. CONCLUSIONS 

The growing attention of the scientific community in 
understanding the behavior of extreme values have led, in 
the recent past, to the development of an extreme value 
theory for dynamical systems. In this framework, it has 
been shown that the statistics of extreme value can be 
linked to the statistics of return in a neighborhood of a 
certain initial conditions by choosing special observables 
that depend on the distance between the iterated trajec- 
tory and its starting point. Until now, rigorous results 
have been obtained assuming the existence of an invari- 
ant measure for the dynamical systems and the fulfill- 
ment of independence requirement on the series of max- 
ima achieved by imposing D' and D 2 mixing conditions, 
or, alternatively, assuming an exponential hitting time 
statistics [22-25]. The parameters of the GEV distri- 
bution obtained choosing as observables the function g^, 
i = 1, 2, 3 defined above depend on the local dimension of 
the attractor D and numerical algorithms to perform sta- 
tistical inference can be set up for mixing systems having 
both absolutely continuous and singular invariant mea- 
sures [27, 28, 39, 40]. Instead, when considering systems 
with regular dynamics, the statistics of the block maxima 
of any observable does not converge to the GEV family 
[18, 19]. 

Taking a complementary point of view, in this paper 
we have studied the statistics of exceedances for the same 
class of observables and derived the limiting distributions 
assuming only the existence of an invariant measure and 
the possibility to define a local dimension D around the 
point £ of interest. To prove that the limiting distribution 
is a GPD we did not use any further conditions. In partic- 
ular no assumptions on the mixing nature of the maxima 
sequence have been made. This means that a GPD lim- 
iting distribution holds for the statistics of exceedancc 
of every kind of dynamical systems and it depends only 
on the threshold value and on the local dimension once 
we choose the observables gi, i = 1, 2, 3. Other functions 



can converge to the limiting behaviour of the GPD family 
if they asymptotically behave like the g^s (compare the 
discussion in [23]). Nonetheless, this requires, analyti- 
cally, to perform separately the limit for the threshold T 
going to g max and that for the radius of the ball going to 
zero. In practical terms, this requires, potentially, much 
stricter selection criteria for the exceedances when finite 
time series are considered. 

We note that, as the parameter £ is inversely propor- 
tional to D, one can expect that each time we analyse 
systems of intermediate or high dimensionality, the dis- 
tributions for gi and g% observables will be virtually in- 
distinguishable from what obtained considering the g\ 
observable: the £ = is in this sense an attracting mem- 
ber of the GPD family. This may also explain, at least 
qualitatively, why the Gumbel (k = for the GEV fam- 
ily) distribution is so efficient in describing the extremes 
of a large variety of natural phenomena [3] . 

The universality of this approach allows to resolve the 
debate on whether there exists or not a general way to ob- 
tain information about extreme values for quasi-periodic 
motions raised in Balakrishnan et al. [18] and Nicolis 
et al. [19]. This is due to the fact that considering several 
of the largest order statistics instead of just the largest 
one we can study orbits where numerous exceedance are 
observed in a given block, as it happens for systems with 
periodic or quasi-periodic behaviors. As an example, one 
may consider a system with multiple commensurable fre- 
quencies: choosing a block length larger or equal to the 
smallest common period, we select always the same value 
for any considered observable. On the other hand for 
mixing maps we find, as expected, an asymptotic equiv- 
alence of the results obtained via the Gnedenko and via 
the Pareto approach. 

The Pareto approach provides a way to reconstruct lo- 
cal properties of invariant measures: once a threshold is 
chosen and a suitable exceedance statistics is recorded, 
we can compute the local dimension for different initial 
conditions taken on the attractor. This is true also in 
the opposite direction: if the knowledge of the exact value 
for the local dimension is available, once we chose a small 
enough radius (threshold), it is possible to compute a pri- 
ori the properties of extremes without doing any further 
computations. In fact, the expression for the parameters 
Eqs. (15)- (17) do not contain any dependence on the 
properties of the dynamics except the local dimension. 

Besides the analytical results, we have proved that 
Pareto approach is easily accessible for numerical investi- 
gations. The algorithm used to perform numerical simu- 
lations is versatile and computationally accessible: unlike 
the GEV algorithm that requires a very high number of 
iterations to obtain unbiased statistics, using the Pareto 
approach we can fix a priori a value for the threshold and 
the number of maxima necessary to construct the statis- 
tics. With the simulations carried out on the standard 
maps, we obtain meaningful results with a much smaller 
statistics with respect to what observed when considering 
the Gnedenko approach. 
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We hope that the present contribution may provide a 
tool that is not only useful for the analysis of extreme 
events itself, but also for characterising the dynamical 
structure of attractors by giving a robust way to com- 
pute the local dimensions, with the new possibility of 
embracing also the case of quasi-periodic motions. 

Future investigations will include the systematic study 
of the impact on the extreme value statistics of adding 
stochastic noise to regular and chaotic deterministic dy- 
namical systems, and the use of the Ruelle response the- 
ory [41, 42] to study the modulation of the statistics of 
extremes due to changes in the internal or external pa- 



rameters of the system, especially in view of potentially 
relevant applications in geophysics such as in the case of 
climate studies [43, 44]. 
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FIG. 1. GPD parameters for the observables gi, i = 1, 2, 3 computed over orbits of the standard map, for various values of the 
constant K. For each value of K, results refer to an esnsemble of 200 randomly chosen initial conditions f . The notation p(gi) 
indicates the parameter p computed using the extreme value statistics of the observable gi. a) £(<7i) VS K, b) £(32) VS K, c) 
£(53) VS K, d) <t(<ji) VS K, e) <t(<?2) VS JiT, f) cr(g3) VS if. Black solid lines: ensemble-average value. Black dotted lines: 
ensemble spread evaluated as one standard deviation of the ensemble. Green lines: theoretical values for regular orbits. Red 
lines: theoretical values for chaotic orbits. 



